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Abstract 


This  report  presents  the  results  of  a  study  undertaken  to  help  build  a  database  and  modeling 
capability  that  will  support  the  design  of  a  laser-based  ignition  system  for  a  gun-fire  simulation 
(GUFS)  platform.  In  this  study,  the  thermal  transport  properties  of  pressed  pellets  of  candidate 
igniter  materials  were  experimentally  investigated  in  the  range  20-100°  C  via  (1)  the  technique 
developed  by  Miller  and  Kotlar  (“Technique  for  Measuring  Thermal  Diffusivity/Conductivity 
of  Small  Thermal  Insulator  Specimens,”  iJeviVw  of  Scientific  Instruments,  vol.  64,  p.  2954, 1 993) 
and  (2)  a  small  modification  to  their  technique  that  allows  it  to  be  used  for  (slightly)  electrically 
conductive  materials.  The  compounds  investigated  in  this  manner  include  black  pow^r,  red 
powder  (with  and  without  0.5-weight-percent  carbon  black),  B/KNO3,  AlA,  a/KClO^, 
Ti/KC104,  an  ARCO  formulation  (with  0.5-weight-percent  carbon  black),  and  an  FeO/Ti/Zr/JNC 
mixture.  In  addition,  the  heat  capacities  of  these  compounds  were  detemmed  over  the  range 
20-100°  C  via  differential  scanning  calorimetry.  The  experimentally  derived  parameters  are 
employed  as  input  to  a  simplified  model  of  the  heat-transfer  process  attending  the  irradiation  of 
a  candidate  with  a  relatively  low-fluence  laser  beam.  The  results  indicate  Aat  the  differences  m 
the  thermal  transport  properties  of  the  materials  studied  will  not  lead  to  significant  differences 
in  the  irradiance  required  to  ignite  them. 
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1.  Introduction 


Laser-based  ignition  is  being  proposed  as  a  means  of  improving  the  safety  of  rounds 
employed  with  a  gun  fire  simulator  (GUFS)  system.  The  system  platform  carries  20  flash- 
and/or  boom-producing  rounds  that  produce  threat  cues  and  target  “hit  and  kill”  recognition 
during  training  exercises.  The  current  rounds  are  initiated  via  electric  matches,  and,  as  with  all 
such  igniters,  they  are  susceptible  to  inadvertent  ignition  by  electrostatic  discharge  (BSD)  and 
electromagnetic  radiation  (ER).  A  number  of  incidents  involving  GUFS  platform  rounds  have 
been  reported,  and  the  risk  that  these  munitions  may  inadvertently  ignite  will  become 
increasingly  acute  as  more  and  more  ER-based  sensors  probe  the  modem  battlefield  and  (by 
extension)  training  exercises.  Laser-based  ignition  systems  are  immune  to  ER  effects,  and  back- 
fitting  the  GUFS  platform  with  such  a  system  would  reduce  the  safety  hazard  associated  with  the 
rounds.  However,  though  there  is  little  doubt  that  a  workable  laser-based  ignition  system  could 
be  built,  the  cost  to  do  so  is  a  concern.  (Current  guidance  is,  that  to  be  tenable,  the  cost/round 
must  be  kept  near  current  levels.)  The  work  reported  here  represents  part  of  an  effort  to  develop 
a  database  and  modeling  capability  that  will  facilitate  the  design  of  a  cost-effective  system. 

Design  options  being  considered  for  the  GUFS  application  include  (1)  having  a  dedicated 
laser  diode  for  each  round  on  the  platform  (20  lasers  total)  or  (2)  having  a  single  laser  diode  (or  a 
solid-state  laser)  coupled  to  each  round  through  an  optical  switch  (one  laser,  one  optical  switch). 
Regardless  of  the  approach  chosen,  the  cost  of  the  system  will  be  driven  by  requirements  for 
laser  power,  optical  component  quality,  and  mechanical  tolerances — all  of  which  depend  on  the 
minimum  irradiance  (intensity/spot  size)  necessary  to  properly  function  a  round.  The  minimum 
irradiance,  in  turn,  depends  on  the  time  available  between  the  arrival  of  a  trigger  signal  and  the 
ignition  event  (i.e.,  ignition  delay)  and  the  mechanism(s)  relevant  to  ignition  on  that  time  scale. 
In  the  case  of  GUFS  platform  rounds,  ignition  delays  of  up  to  0.5  s  are  permissible.  This  time 
scale  allows  relatively  low  laser  fluences  to  be  employed,  and  the  fluence  will  act  simply  as  a 
heat  (flux)  source  that  raises  the  temperature  of  the  igniter  material  until  self-propagating 
chemical  reactions  ensue. 
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A  simplified,  one-dimensional  (1-D)  governing  equation  for  the  spatial  (x)  and  temporal  (t) 
dependence  of  the  temperature,  T(x,t),  in  a  material  under  such  flux  conditions  is  [1], 

K  =  +  Qiexp(.-E/RT).  (1) 

dt  dx^  pcp  dx  Cp 

The  first  term  on  the  right-hand  side  of  this  equation  is  associated  with  the  diffusion  of  heat  due 
to  temperature  gradients,  the  dynamics  of  which  depend  on  the  material  s  thermal  diffusivity  (OJ). 
(The  thermal  conductivity  [A]  of  a  material  is  related  to  its  thermal  diffusivity  per, 


A 


a  = - , 


PCp 


(2) 


where  p  is  the  material’s  density  and  Cp  is  its  heat  capacity  at  constant  pressure.)  The  second 
term  on  the  right-hand  side  of  equation  (1)  is  associated  with  the  absorption  of  radiant  energy  by 
the  material,  the  dynamics  of  which  depend  on  the  incident  flux  (fp)  and  the  material  s  absorption 
coefficient  (n).  The  last  term  is  associated  with  chemical  energy  production,  with  the  given 
function  being  the  Arrhenius  form  of  the  heat  release  rate.  The  parameters  of  this  function  are 
the  activation  energy  (E),  a  first-order  pre-exponential  factor  (z),  and  the  heat  of  reaction  (Q).  R 
is  a  constant. 

Unfortunately,  even  this  simplified  model  requires  a  plethora  of  parameters  that  are  not 
generally  available,  making  it  difficult  to  identify  systems  with  potentially  low-fluence 
requirements.  The  work  presented  in  this  report  is  concerned  with  the  thermal  transport 
properties  of  various  igniter  materials  and  the  influence  thermal  diffusion  will  have  in 
determining  the  minimum  fluence  required  to  initiate  these  materials.  Presented  are  the  results 
of  experiments  conducted  to  establish  the  thermal  transport  properties  of  nine  pyrotechnic 
compounds  over  the  range  20-100°  C.  Most  measurements  were  made  with  a  slightly  modified 
version  of  the  experimental  technique  developed  by  Miller  and  Kotlar  (MK)  [2],  the  modified 
experiment  yielding  a  measurement  for  the  parameter  (XspsCpsf^.  hi  the  interest  of  establishing 
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A  and  a  from  these  measurements,  heat  capacities  for  all  of  the  compounds  were  determined  for 
the  range  20-100°  C  via  differential  scaiming  calorimetry.  To  simulate  the  temperature  response 
of  the  candidates  to  relatively  low-level  laser  fluences,  the  experimental  values  were  employed 
as  parameters  in  a  1-D  heat-transfer  model  that  assumes  n  =  <»  and  Q  =  0,  and  the  irradiance 
required  to  raise  the  compound’s  surface  to  its  autoignition  temperature  within  0.5  s  is 
calculated.  The  results  indicate  that  the  differences  in  the  thermal  transport  properties  of  the 
candidate  materials  (alone)  will  not  lead  to  significant  differences  in  the  irradiance  required  to 
ignite  them. 

2.  Description  of  Igniter  Materials 

Table  1  shows  the  chemical  compositions  and  some  of  the  physical  properties  of  the  igniter 
materials  identified  by  the  U.S.  Army  Armament  Research,  Development,  and  Engineering 
Center  (ARDEC)  as  candidates  for  the  GUFS  application.  Except  for  red  powder,  all  are 
“standard”  pyrotechnics  for  which  some  degree  of  characterization  exists  in  the  open  literature. 
Red  powder  is  a  patented  ARDEC  formulation  that  is  being  developed  in  conjunction  with 
Thiokol.  Lots  M292  and  M296  were  doped  with  a  small  amount  of  carbon  black.  It  is 
considered  that  such  an  additive  may  reduce  the  fluence  required  to  ignite  a  material  by 
increasing  its  optical  density.  However,  it  remains  to  be  determined  whether  the  dopant 
produces  other  effects,  such  as  increasing  the  thermal  conductivity,  that  negate  this  potential 
benefit. 

All  compounds  employed  in  the  experiments  were  prepared  by  ARDEC  and  provided  to  us 
as  cylindrical  pellets  (0.65-cm  diameter  x  0.65-cm  length).  The  pellets  were  formed  by  pressing 
the  mixtures  to  152  MPa  (22  kpsi).  The  pellets  were  friable,  particularly  those  with  high  metal 
or  metal  oxide  content,  and  care  needed  to  be  exercised  in  handling  them.  For  most  of  the 
measurements  reported  here,  the  pellets  were  simply  used  as  supplied,  fri  cases  where  thinner 
specimens  were  employed  (i.e.,  for  the  straight  adaptation  of  the  MK  technique  [2]),  the  pellets 
were  wafered  by  inserting  them  into  an  acrylic  tube  and  cutting  the  resulting  assembly  with  a 
low-speed  saw  with  a  diamond-coated  blade.  In  all  of  these  cases,  ethyl  alcohol  was  employed 
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Table  1.  Pyrotechnic  Compositions  and  Some  of  Their  Properties 


Lot  No.“ 

Material*’ 

(weight-percent  components) 

Pellet  Loading 
(MPa) 

Autoignition 

Temperature 

(“C) 

M279 

Red  Powder 

(PhPhK-1 5/KNO3-75/S-IO) 

152 

1.86 

450 

M291 

Zr/KClOVVAAR 

(Zr-50/KC104-50/VAAR-l) 

152 

2.86 

41  r 

M292 

ARCO  Powder  w/  0.5%  Carbon 

Black 

152 

1.74 

Unknown 

M293 

Black  Powder 
(C-I5/KNO3-75/S-IO) 

152 

1.74 

464^* 

200'* 

M294 

B/KNO3 

(B-25/KNO3-75/V  AAR-1 ) 

152 

1.59 

431 

M295 

Ti/KClOVVAAR 

(Ti-50/KC104-50/VAAR-l) 

152 

2.34 

Unknown 

M296 

Red  Powder  w/  0.5%  Carbon  Black 
(PhPhK-15/KNO3-75/S-10) 

152 

1.85 

450 

M297 

FeO/Ti/Zr/NC 

(FeO-50/Ti-32.5/Zr-17/NC-5) 

152 

3.03 

456 

SI163 

AlA 

(Zr-65/Fe2O3-25/DE-10/VAAR-l) 

152 

3.18 

427 

®  ARDEC  designation. 

PhPhK  =  potassium  salt  of  phenolphthalein,  DE  =  diatomaceous  earth,  and  VAAR  =  vinyl  alcohol  acetate  resin. 
Value  is  the  material’s  explosive  temperature,  the  autoignition  temperature  is  unknown. 

Two  different  autoignition  values  are  reported  for  black  powder. 


as  a  cutting  fluid.  In  some  cases,  the  specimen’s  produced  by  wafering  were  sanded  with  IS-pm 
grit  paper  to  improve  their  surface  planarity  and  finish.  If  wafered,  the  specimens  were  stored  in 
an  evacuated  dessicator  for  at  least  2  days  prior  to  using  them  in  an  experiment. 


3.  Experimental  Considerations 


Except  for  minor  variations,  the  experimental  apparatus  used  in  this  study  was  the  same  as 
that  developed  and  employed  by  MK  [2].  MK  have  discussed  the  difficulty  of  measuring  the 
thermal  conductivity  and  thermal  diffusivity  of  small  solid  samples,  and  their  technique 
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addresses  t5^ical  concerns  such  as  the  generation  of  a  well-defined  heat  flux,  the  match  between 
experimental  and  model  boundary  conditions,  and  instrumentation  errors.  A  schematic  diagram 
of  their  experiment  is  shown  in  Figure  1.  The  experiment  involves  resistively  heating  a 

4.4- cm-long,  1.2-cm-wide,  0.0005-cm-thick  Constantan  foil  (Hamilton  Precision  Metals)  via  a 
voltage-regulated  power  supply  (Hewlett-Packard,  Model  6024A),  the  circuit  producing  a  step 
increase  in  current  (and  heat  flux)  when  completed  via  a  mercury-wetted  relay.  The  heat  flux 
sends  thermal  waves  into  the  adjacent  test  sections,  and  the  temperature  response  at  a  known 
distance  from  the  foil  is  monitored  with  a  bare  Chromel/Alumel  thermocouple  (Omega,  Style  I 
or  Style  II).  The  output  of  the  thermocouple  is  preamplified  via  a  differential  amplifier  (Stanford 
Research  Systems,  Model  SR560)  and  recorded  by  a  digital  oscilloscope  (Nicolet,  Model 
4094/4570)  over  a  4-  or  8-s  period.  One-tenth  of  the  period  is  recorded  prior  to  heat  flux  onset, 
providing  an  unbiased  window  on  the  test  fixture  temperature  (To).  A  data  point  was  acquired 
every  0.002  s  for  the  4-s  duration  experiments  and  every  0.004  s  for  the  8-s  duration 
experiments.  To  determine  the  heat  flux  (fo)  generated  by  the  foil,  a  nominal  value  for  the 
current  in  the  foil  (/«,)  was  obtained  from  the  readout  of  a  multimeter  (Fluke,  Model  77)  and  the 
(fvf)  heat  generation  was  calculated  using  a  previously  measured  value  of  the  foil’s  resistance 
(r/  =  0.34  Q).  The  power  supply  was  typically  set  to  produce  a  current  level  in  the  range 

1.5- 2.5  A,  with  the  resultant  heat  flux  being  in  the  range  0.14-0.40  W/cm^.  The  temperature  of 
the  test  fixture  is  controlled  by  circulating  a  50/50  ethylene  glycol-water  mixture  from  a 
temperature-controlled  bath  (Neslab,  Model  ULT-80DD)  through  copper  tubing  soldered  to  the 
copper  block  housing  the  test  fixture. 


Material  with  unknown 
thermal  transport  properties 


Thermocouple  _ 

location 

Material  with  known 
thermal  transport  properties 


Thermal  guard 


SAMPLE  SIDE(x>0) 


Heat’source  plane  (x  *  0) 
(Constantan  foil) 


BASE  SIDE(x<0) 


Figure  1.  Schematic  Diagram  of  Experimental  Configuration  Developed  by  MK  [2]. 
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MK  have  shown  that  their  experiment  approaches  a  mathematically  idealized  system  in 
which  a  step  function  heat  flux  emanates  from  a  plane  between  two  semi-infinite  media.  The 
governing  equation  for  this  idealization  is 


dT 

—  =  a-^, 

dt  dx 


(3a) 


with  boundary  conditions 


T(x,0)  =  Ty,-co<  X  <+oo, 


(3b) 


dT(O-,0 

dx 


=  /  (0  ,t);  x<0. 


(3c) 


^  =  /^(o",t);x>0,and 

OX 


(3d) 


T{o\t)  =  no-,t), 


(3e) 


where  x  is  the  distance  from  the  source  plane  and  the  subscripts  s  and  b  identify  properties 
belonging  to  the  sample  and  base,  respectively.  Conservation  of  energy  requires  that  the  sum  of 
the  heat  fluxes  going  into  the  sample,  /^(0^^i  and  the  base, /'(O’,?)  ,  be  equal  to  the  total  flux 
emanating  from  the  source  plane,  /».  The  solution  to  this  system  of  equations  [i.e.,  the 
temperature  response  of  a  point  in  the  sample  (x  >  0)]  is  given  by  [3] 


The  last  term  on  the  right-hand  side  of  equation  (4)  is  the  complementary  error  function  integral 
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(5) 


To  derive  ctj  and  equation  (4)  is  fit  to  the  experimentally  measured  temperature  response 
via  a  nonlinear  least-squares-fitting  routine.  (The  computation  of  the  complex  error  function 
integral  was  accomplished  by  calculating  the  complex  error  function  for  the  argument  via  the 
rational  approximation  suggested  by  Hui  et  al.  [4]  and  making  an  appropriate  transformation.)  In 
performing  a  fit,  all  model  parameters  except  (Xs,  Aj,  and  To  are  set  to  measured  values  and  the 
fitting  routine  allowed  to  vary  ctj,  A®,  and  To  to  minimize  the  root  mean  square  of  the  residual. 
Being  mathematically  independent  parameters,  both  cCs  und  A^  are  determined  directly  from  this 
experiment. 

It  was  found,  however,  that  the  MK  technique  could  not  be  used  to  characterize  the  thermal 
transport  properties  of  the  FeO/Ti/Zr/NC  compound  [5].  During  experiments  with  this 
compound,  the  thermocouple  reading  dropped  precipitously  at  the  instant  the  electrical  heating 
pulse  started.  This  suggested  that  current  was  leaking  to  the  thermocouple,  compromising  not 
only  the  thermocouple  output,  but  our  ability  to  characterize  the  heat  flux.  This  observation 
prompted  us  to  consider  monitoring  the  temperature  response  in  the  base.  (See  Figure  2.)  Not 
only  does  this  configuration  resolve  the  problem  posed  by  (slightly)  electrically  conductive 
materials  such  as  FeO/Zr/Ti/NC,  it  offers  operational  advantages.  It  precludes  the  need  to 
prepare  thin  samples,  and  it  facilitates  the  testing  of  a  series  of  samples  because,  once  in  place, 
the  fragile  thermocouple  does  not  have  to  be  handled  again. 

However,  the  benefits  accrued  through  this  approach  come  with  a  significant  price — ^namely, 
that  the  function  describing  the  temperature  response  becomes 


T(x,t)  =  To 


(6a) 


or 


7 


SAMPLE  SIDE(x>0) 


Material  with  unknown 
thermal  transport  properties 


Thermocouple  _ _ 

location 

Material  with  known 
thermal  transport  properties 

Thermal  guard 


Heat  source  plane  (x  0) 
(Constantan  foil) 


BASE  SIDE(x<0) 


Figure  2.  Schematic  Diagram  of  Experimental  Setup  Employed  for  This  Study. 


T{x,t)=T„ 


^^Pb^pb 


ierfc 


(6b) 


and  in  fitting  this  function  to  the  experimental  data,  a,  and  Is  are  no  longer  mathematically 

1/2 

independent  parameters.  Instead,  the  parameter  derivable  from  the  experiment  is  (XspsCps)  . 
Thus  a  specimen’s  p  and  Cp  must  be  known  to  get  absolute  values  for  Xs  and  While  the 
density  is  a  relatively  easy  property  to  measure,  heat-capacity  measurements  are  a  nontrivial 
extra  step.  Nevertheless,  the  Figure  2  configuration  was  used  for  most  of  the  results  reported 
here.  To  determine  (XspsCpf^  from  experiments  with  this  configuration,  nonlinear  least-squares 


fits  of  equation  (6b)  to  the  experimental  data  were  performed,  all  model  parameters  except 
(XsPsCps)^'^  and  To  being  set  to  measured  values,  and  the  program  allowed  to  vary  (XspsCps)^’'^  and 


To  to  optimize  the  fit  of  the  function  to  the  data. 


Another  change  we  made  to  the  basic  technique  developed  by  MK  involved  refining  the 
temperature  response  model  to  account  for  the  time  constant  of  the  thermocouple.  This  change 
was  necessitated  by  the  use  of  Style  I  thermocouples  instead  of  Style  H  thermocouples  to 
monitor  the  temperature  at  one  juncture  in  this  study.  Though  made  from  the  same  thickness  of 
foil.  Style  I  and  Style  n  thermocouples  have  different  geometries,  and  the  time  constant  of  Style 
I  thermocouples  (tc  =10-20  ms)  is  longer  than  that  of  Style  n  thermocouples  (4  =  2-5  ms). 
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(The  seller  of  the  thermocouples  defines  the  time  constant  simply  as  the  time  required  for  the 
thermocouple  reading  to  transition  63.2%  of  an  instantaneous  temperature  change.)  The 
difference  in  time  constants  proved  to  be  significant,  with  <x  and  A  values  derived  from 
experiments  with  Style  I  and  Style  n  thermocouples  being  different  if  the  time  constant  was 
neglected.  This  is  in  contrast  to  the  case  examined  by  MK,  who  found  that  the  results  obtained 
using  5-  and  12-pm-thick  (Style  II)  thermocouples  were  the  same  and  concluded  that  the  time  lag 
of  the  thermocouple  response  could  be  neglected.  To  incorporate  the  effect  of  the  time  constant 
into  the  data  reduction  model,  it  was  assumed  that  the  response  of  the  thermocouple  iR(t))  to  a 
unit  step  increase  in  temperature  was  associated  with  the  form 

Rit)  =  1  -  exp(-t/tc)  (7) 

and  EhihameTs  formula  [6]  was  employed  to  express  the  response  of  the  measurement 
system— the  oscilloscope  voltage  reading,  u(,t)— to  the  driving  function  T(x,t)  [i.e.,  equation  (4) 
or  (6b)]: 


u(jt)  =  t  R’('zyr(.x,t-x)dx.  (8) 

Jo 

This  equation  was  incorporated  into  the  system  model  by  calculating  R\i)T{x,t  -  t)  at  0.001 -ms 
intervals  for  0  <  x  <  t  and  using  Simpson’s  rule  to  evaluate  the  integral.  Corrected  in  this 
manner,  the  results  obtained  with  Style  I  and  Style  II  thermocouples  were  brought  into 
agreement. 

It  is  also  worth  noting  that  fused  silica  was  employed  as  the  base  material  in  most  of  the 
experiments.  MK  had  employed  acrylic  as  the  base  material,  and  the  initial  room-temperature 
results  of  this  study  were  obtained  with  an  acrylic  base.  However,  acrylic’s  glass  transition 
temperature  is  below  100°  C,  a  value  that  we  were  interested  in  reaching.  The  decision  to 
employ  fused  silica  was  based  primarily  on  its  availability  and  having  had  previous  experience 
with  it  in  this  experiment.  But,  it  is  not  an  optimal  choice.  As  noted  by  MK,  the  best 
measurements  will  be  obtained  when  the  thermal  transport  properties  of  the  base  and  sample  are 
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similar.  [It  is  suggested  that  {XspsCps)^'^  be  used  in  matching  materials  for  this  application.]  The 
literature  value  of  iXsPsCp^^'^  for  fused  silica  is  about  0.14  J/cm^-K-s'^^,  while,  for  the  compounds 
studied  here,  this  parameter  was  found  to  fall  in  the  range  0.05-0.10  J/cm  -K-s  . 

The  values  of  a  and  A  for  the  base  material,  which  are  needed  as  input  to  the  model  of  the 
heat-transfer  process,  were  established  experimentally  by  conducting  tests  in  which  the  sample 
side  was  occupied  by  a  specimen  from  the  same  lot  as  that  used  to  construct  the  base.  This 
configuration  provides  a  direct,  simultaneous  measurement  of  a  and  A  without  any  a  prion 
knowledge  of  these  properties.  Values  for  the  fused  silica  base  were  obtained  for  the  range 
20-100°  C  and  the  linear  least-squares  fit  of  this  data — X  =  12.055  -  8.4  x  10  ^TmW/cm-K,  a  = 
7.6957  -  2.14  X  lO’^r  cm^/s,  where  T  is  in  degrees  Celsius — ^was  incorporated  directly  into  the 
data-reduction  model.  These  values  are  similar  to  those  reported  by  MK  [2]. 


Toward  establishing  a  and  A  for  the  pyrotechnic  compounds  from  {XsP^p^^'^  values, 
heat-capacity  measurements  were  made  with  a  Mettler  TA3000  differential  scanning  calorimeter 
(DSC).  DSC  measurements  of  heat  capacity  are  based  on  comparing  the  electrical  power  needed 
to  ramp  (at  a  constant  rate)  the  temperature  of  an  empty  (aluminum)  pan  and  the  same  pan  with  a 
known  mass  of  sample.  For  the  results  reported  here,  a  sample  of  10  to  15  mg  was  shaved  from 
a  pellet,  introduced  to  the  DSC,  cooled  to  -30°  C,  then  ramped  at  10°  C/min  up  to  1 10°  C.  After 
this  measurement,  the  sample  was  “dried”  by  heating  it  at  100°  C  for  approximately  1  hr  with 
5-10  ml/min  of  argon  flowing  though  the  chamber.  The  sample  was  then  cooled  back  down, 
reweighed,  and  the  measurement  was  repeated.  Except  where  noted,  the  values  reported  in 
section  4  are  derived  from  the  second  measurement  on  each  sample.  The  values  include  a 
correction  based  on  a  temperature-dependent  scaling  factor  derived  from  tests  with  sapphire 
specimens.  (The  scaling  factor  adusts  the  DSC  results  for  the  sapphire  specimens  such  that 
literature  values  are  obtained.) 

For  the  purpose  of  deriving  A  and  (X  from  {X^PsCps)  » the  results  for  two  different  samples 
were  averaged  and  a  polynomial  of  the  form 

Cp(T)  =  ao  +  aiT  +  a2T^  +  +  ...  (9) 
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fit  to  the  average.  The  order  of  the  polynomial  employed  to  fit  the  data  was  based  on  our 
judgement  as  to  a  polynomial’s  ability  to  reproduce  the  data  over  the  temperature  range  of 
interest  and  on  using  the  lowest  order  polynomial  possible. 


It  is  also  noted  that  the  DSC  measurements  (should  have)  provided  Cp  values  for  the  range 
-10-100°  C.  However,  in  all  of  the  compounds  studied,  the  measured  Cp  values  dropped 
relatively  rapidly  with  temperature  for  temperatures  below  20°  C.  This  trend  was  at  odds  with 
expectations  for  compounds  whose  Cp  values  could  be  estimated  based  on  literature  values  for 
their  constituents,  prompting  us  to  develop  an  independent  experimental  check  of  the  DSC 
results.  The  check  involved  calculating  Cp  for  the  fused  silica  employed  as  the  base  in  the  test 
fixture  from  equation  (2)  using  (1)  values  of  a  and  X  found  via  the  MK  technique  and  (2)  a 
separate  measurement  of  density.  As  shown  in  Figure  3,  the  values  derived  from  the  two 
techniques  are  in  good  agreement  for  the  range  20-100°  C.  They  are,  however,  significantly 
below  values  recommended  in  the  literature  [7].  Thus,  though  we  have  some  misgivings  about 
the  Cp  results,  the  fused  silica  results  and  comparison  with  estimates  for  some  of  pyrotechnic 
compounds  based  on  literature  values  for  their  constitutents  support  their  validity  in  the  range 
20-100°  C.  Also,  since  this  range  covers  most  of  the  thermal  transport  data  collected,  no  effort 
was  undertaken  to  reconcile  the  DSC  Cp  results  and  expectations  for  lower  temperatures  and  the 
lower  temperature  results  are  not  reported  here. 


- Ref.  7 

♦  MK  meas 

—  DSC  (A) 

—  DSC  (B) 


temperature  (C) 

Figure  3.  Comparison  of  Experimentally  Measured  Values  of  the  Heat  Capacity  of  Fused 
Silica  and  the  Recommendation  of  Reference  [7]. 
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4.  Experimental  Results 


Figures  4-12  show  the  experimentally  determined  heat  capacities  for  the  nine  candidate 
compounds  provided  to  us  by  ARDEC.  In  cases  where  the  heat  capacities  of  the  constituents 
could  be  found  in  the  literature  [7],  an  estimate  for  the  compound  is  provided  for  comparison. 
(The  mass  fractions  of  the  constituents  needed  for  these  estimates  are  given  in  Table  1.)  In  the 
case  of  the  compounds  containing  red  powder  (M279  and  M296),  the  Cp  values  obtained  from  the 
first  transient  to  110°  C  are  higher  than  those  measured  during  the  second  transient.  It  was  also 
found  in  these  two  cases  that  the  mass  of  the  sample  was  reduced  on  the  order  of  1%  by  the  first 
measurement  and  drying  period.  (The  mass  loss  produced  by  the  second  transient  was 
negligible.)  A  comparison  of  the  results  (Figures  6  and  7)  suggests  that  a  gas-producing 
“reaction,”  perhaps  simply  evaporation,  onsets  at  about  50°  C  in  these  compounds. 


temperature  (C) 


Figure  4.  DSC  Measurements  of  the  Heat 
Capacity  of  Black  Powder  for 
Temperatures  From  20  to  100°  C. 


Figure  5.  DSC  Measurements  of  the  Heat 
Capacity  of  B/KNO3  for 
Temperatures  From  20  to  100°  C 
and  an  Estimate  Based  on 
Literature  Values  for  the 
Constituents. 
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temperature  (C) 

Figure  6.  DSC  Measurements  of  the  Heat 
Capacity  of  Red  Powder  for 
Temperatures  From  20  to 
100°  C.  Comparison  of  Results 
Obtained  Before  and  After 
Drying. 


temperature  (C) 

Figure  8.  DSC  Measurements  of  the  Heat 
Capacity  of  Zr/KC104  for 
Temperatures  From  20  to  100°  C 
and  an  Estimate  Based  on 
Literature  Values  for  the 
Constituents. 


temperature  (C) 

Figure  7.  DSC  Measurements  of  the  Heat 
Capacity  of  Red  Powder  With 
0.5-Weight-Percent  Added 
Carbon  Black  for  Temperatures 
From  20  to  100°  C.  Comparison 
of  Results  Obtained  Before  and 
After  Drying. 


20  30  40  50  60  70  80  90  100 


temperature  (C) 

Figure  9.  DSC  Measurements  of  the  Heat 
Capacity  of  Ti/KC104  for 
Temperatures  From  20  to  100°  C 
and  an  Estimate  Based  on 
Literature  Values  for  the 
Constituents. 


13 


heat  capacity  (J/g-C) 


20  30  40  50  60  70  80  90  100 


temperature  (C) 

Figure  10.  DSC  Measurements  of  the  Heat 
Capacity  of  FeO/Ti/Zr/NC  for 
Temperatures  From 
20  to  100°  C  and  an  Estimate 
Based  on  Literature  Values  for 
the  Constituents. 


temperature  (C) 


Figure  11.  DSC  Measurements  of  the  Heat 
Capacity  of  AlA  for 
Temperatures  From 
20  to  100°  C. 


20  30  40  50  60  70  80  90  100 


temperature  (C) 

Figure  12.  DSC  Measurements  of  the  Heat  Capacity  of  an  ARCO  Formulation  With 

0.5- Weight-Percent  Added  Carbon  Black  for  Temperatures  From  20  to  100°  C. 
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Figure  13  shows  the  values  of  (IsPsCps)^^^  for  six  of  the  nine  candidate  compounds  as 
determined  from  the  modified  version  of  the  MK  technique.  Values  are  not  provided  for  the  two 
compounds  with  red  powder  (M279  and  M296)  or  A1 A  (SI-163).  Experiments  were  conducted 
with  specimens  containing  red  powder,  but  the  results  of  the  heat-capacity  measurements,  which 
show  that  their  properties  depend  on  their  thermal  history,  argue  against  reducing  the  data  and 
reporting  the  results.  The  AlA  pellets  proved  to  be  very  friable,  and  none  of  the  specimens 
ended  up  being  of  sufficient  quality  to  conduct  the  thermal  transport  property  measurements. 
(The  better  pellets  that  we  received  were  consumed  in  attempts  to  wafer  them.)  The  reported 
values  are  observed  to  fall  within  a  fairly  narrow  range  that  is  bounded  above  by  those  for 
FeO/Ti/Zr/NC  and  below  by  those  for  Zr/KC104.  The  values  are  similar  to  that  of  acrylic  (0.053 
J/cm^-K-s^'^)  [2].  In  aU  cases,  (XspsCpsf^  is  observed  to  increase  slightly  with  temperature  over 
the  range  20-100°  C. 

Figures  14-19  show  the  temperature-dependent  values  of  a  and  A  for  the  six  compounds  as 
derived  from  the  results  with  the  modified  (Figure  2)  experimental  configuration  and  the  DSC 
heat-capacity  measurements.  A  surprising  result  was  the  difference  in  property  values  obtained 
in  this  manner  and  values  obtained  earlier  with  the  basic  (Figure  1)  configuration  [5].  (See 
Table  2.)  It  is  assumed  that  the  discrepancies  are  associated  with  the  fact  that  the  compounds  are 
slightly  electrically  conductive  and  that  current  leakage  from  the  foil  to  the  thermocouple 
impacts  the  temperature  reading.  Such  an  effect  was  obvious  in  the  case  of  FeO/Ti/Zr/NC.  It  is 
presumed  that  the  effect  is  just  subtler  for  these  other  compounds. 

5.  Influence  of  Thermal  Transport  Properties  on 
Irradiance  Requirements 

Equation  (4)  describes  the  temperature  transient  that  occurs  during  low-fluence  irradiation  of 
a  pyrotechnic  pellet  in  contact  with  a  flat,  transparent  window  if  (1)  all  of  the  incident  radiation 
is  absorbed  at  the  surface  of  the  pellet  (i.e.,  n  =  «)  and  (2)  energy  production  (or  loss)  due  to 
chemical  reactions  or  phase  changes  is  negligible  {Q  =  0).  Further  noting  that  the 
complementary  error  function  integral  decreases  as  its  argument  increases,  and  thus  that  the 
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B/KN03 
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Formulation  With  0.5-Weight-Percent  Added  Carbon  Black. 


Figure  14.  Thermal  Diffusivity  and 
Thermal  Conductivity 
Measurements  for  Black 
Powder. 


temperature  (C) 


Figure  16.  Thermal  Diffiisivity  and 
Thermal  Conductivity 
Measurements  for  Zr/KC104. 


Figure  15.  Thermal  Difihisivity  and 
Thermal  Conductivity 
Measurements  for  B/KNO3. 


temperature  (C) 


Figure  17.  Thermal  Di^sivity  and 
Thermal  Conductivity 
Measurements  for  Ti/KC104. 
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Figure  18.  Thermal  Diffusivity  and 
Thermal  Conductivity 
Measurements  for  an  ARCO 
Formulation  With 
5-Weight-Percent  Added 
Carbon  Black. 


Figure  19.  Thermal  Diffusivity  and 
Thermal  Conductivity 
Measurements  FeO/Ti/Zr/NC. 


Table  2.  Comparison  of  Thermal  Transport  Property  Measurements  at  23.5°  C 


Material 

Source 

A 

(mW/cm-C) 

Cp 

(J/g-°C) 

Black  Powder 

McQuaid  and  Beyer  [5] 

3.3 

2.5 

0.8 

— 

This  work 

2.0 

1.0 

1.1 

B/KNO3 

McQuaid  and  Beyer  [5] 

3.6 

3.0 

0.8 

— 

This  work 

1.2 

0.7 

1.1 

Zr/KC104 

McQuaid  and  Beyer  [5] 

4.8 

3.4 

0.5 

— 

This  work 

1.9 

1.2 

0.6 

Acrylic 

Miller  and  Kotlar  [2] 

1.8 

1.1 

1.3 

maximum  temperature  will  be  at  the  surface  (x  =  0)  for  all  t,  the  irradiance  necessary  to  raise  the 
temperature  of  the  surface  to  the  material’s  autoignition  temperature  under  this  set  of 
assumptions  may  be  calculated  by  setting  x  equal  to  zero  and  rearranging  equation  (4)  to  give 
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(10) 


where  AT  is  the  difference  between  the  autoignition  temperature  and  the  ambient  temperature. 
Expressed  in  this  manner,  the  desirability  of  having  materials  with  low  (product)  values  of  A,  p, 
and  Cp  is  apparent.  Moreover,  this  expression  shows  that  the  results  presented  in  Figure  13  can 
be  used  directly  in  estimating  flux  requirements  under  the  given  set  of  assumptions. 

To  provide  some  appreciation  of  the  ramifications  of  a  material’s  thermal  transport  properties 
on  the  fluence  required  to  ignite  it,  we  assume  that  the  results  in  Figure  1 3  can  be  extrapolated  to 
the  autoignition  temperature  for  each  compound  and  calculate  the  fluence  needed  to  raise  the 
temperature  of  the  surface  of  each  material  from  25°  C  to  its  autoignition  temperature  in  0.5  s. 
(Again,  0.5  s  is  the  maximum  ignition  time  delay  permitted  for  the  GUFS  application.)  The 
results  for  three  different  cases  are  presented  in  Table  3.  In  Case  A,  it  is  assumed  that  the 
bounding  window  is  a  perfect  insulator  (A*  =  0).  In  Case  B,  the  bounding  window  is  given  the 
properties  of  acrylic.  And,  in  Case  C,  the  bounding  window  is  given  the  properties  of  fused 
silica.  (In  each  case,  the  thermal  transport  properties  of  the  bounding  window  are  assumed  to  be 
temperature  independent.)  For  reference,  130  W/cm^  is  equal  to  10  mW  at  the  distal  end  of  an 
optical  fiber  with  a  lOO-pm-diameter  core — a  value  that  would  make  laser-based  ignition  a 
viable  option  for  the  GUFS  system.  Except  for  the  value  for  black  powder  computed  assuming 
an  autoignition  temperature  of  200°  C,  the  differences  between  compounds  for  a  given  case  are 
not  considered  significant.  The  difference  between  cases  are  notable,  with  the  minimum  fluence 
for  Case  B  being  approximately  50%  higher  than  for  Case  A  and  the  minimum  fluence  for 
Case  C  being  approximately  twice  as  high  as  Case  A. 

6.  Summary 

It  has  been  recognized  that  the  thermal  transport  properties  of  the  igniter  material  in  a 
laser-based  ignition  system  will  influence  the  minimum  irradiance  required  to  function  a  round 
[8].  However,  the  difficulty  of  obtaining  reliable  thermal  transport  property  data, 
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Table  3.  Minimum  Fluence  Required  to  Reach  Autoignition  Temperature  Within  0.5  s 


Lot  No. 

Material 

70 

100 

140 

M293 

Black  Powder 

20 

50 

90 

M294 

B/KNO3 

70 

100 

140 

M291 

Zr/KClOVVAAR* 

50 

80 

120 

M295 

Ti/KC104*’ 

80 

110 

150 

M297 

FeO/Ti/Zr/NC 

60 

90 

130 

M292 

ARCO  w/  carbon  black*’ 

50 

80 

120 

^  Values  for  this  material  are  based  on  its  explosive  temperature. 

^  The  autoignition  temperature  is  unknown.  A  value  of  450®  C  has  been  assumed  for  the  calculations. 


particularly  for  energetic  materials,  is  a  hurdle  to  quantifying  the  impact  of  this  aspect  of  the 
energy  transport  process.  The  results  of  this  study  provide  values  for  the  thermal  transport 
properties  of  six  pyrotechnic  compounds  being  considered  for  use  in  back-fitting  a  GUFS  system 
with  a  laser-based  ignition  system.  The  results  were  obtained  with  an  experimental  technique 
developed  specifically  for  the  difficult  case  of  acquiring  such  information  from  small, 
electrically  conductive  specimens  with  relatively  low  diffiisivity  values.  In  addition  to  this 
primary  objective,  DSC  results  for  compounds  containing  red  powder  indicate  that  the  thermal 
stability  of  this  pyrotechnic  may  be  an  issue. 

To  begin  to  quantify  the  impact  of  thermal  transport  properties  on  irradiance  requirements, 
the  experimental  results  were  employed  as  input  to  a  model  that  assumes  all  the  laser  energy  is 
absorbed  at  the  surface,  no  phase  changes  or  chemical  reactions  occur,  and  that  thermal  transport 
property  measurements  for  each  compound  can  be  extrapolated  to  their  autoignition  temperature. 
While  this  model  is  a  great  oversimplification  of  the  actual  process,  and  thus  of  limited  value  for 
establishing  minimum  irradiance  requirements  for  such  systems,  its  predictions  are  consistent 
with  expectations  based  on  previous  studies.  Moreover,  the  results  obtained  for  the  different 
materials  are  instructive.  In  particular,  they  indicate  that  the  differences  in  thermal  properties  of 
the  candidates  studied  thus  far  will  not  be  a  significant  factor  from  the  standpoint  of  reducing 
fluence  requirements.  On  the  other  hand,  the  packaging  of  the  igniter  material  will  be  an 
important  consideration. 
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List  of  Symbols 


Nomenclature 

a  thermal  diffusivity  (square  centimeters/second) 

A  thermal  conductivity  (watts/centimeter»degrees  Celsius) 

p  density  (grams/cubic  centimeter) 

A  area  of  the  foil  (square  centimeters) 

Cp  heat  capacity  (joules/gram»degrees  Celsius) 

E  activation  energy  in  the  Arrhenius  heat  release  expression  (joules/mole) 
fo  heat  flux  (watts/square  centimeter) 

/■(0,r)  heat  flux  into  base  (watts/square  centimeter) 

f\Q,t)  heat  flux  into  sample  (watts/square  centimeter) 

lo  nominal  current  of  electrical  pulse  employed  to  heat  foil  (amperes) 

I  length  (centimeters) 

n  absorption  coefficient 

Q  heat  of  reaction 

r/  resistance  of  foil  (ohms) 

R  gas  constant  (joules/degrees  Celsius«mole) 

R(t)  temporal  response  of  passive  system  to  a  unit  step  forcing  function 
t  time  (seconds) 

To  ambient  temperature  (degrees  Celsius) 

T(x,i)  actual  temperature  (degrees  Celsius) 
u(t)  measured  temperature  (degrees  Celsius) 

X  distance  from  heat  source  (centimeters) 

z  first  order  pre-exponential  factor  in  the  Arrhenius  heat  release  expression 

Subscripts 

b  base  material  property 

s  specimen  property 
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